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ABSTRACT 



A method is described for determining analytically or graphically, from appro- 
priate time-temperature freezing and melting curves obtained on hydrocarbons, 
the freezing point of a given substance, and within certain wide limits, the freezing 
point of that substance for zero impurity (liquid-soluble, solid-insoluble). 
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I. INTRODUCTION 

In a recent paper [1] 4 dealing with time-temperature freezing and 
melting curves obtained on l^drocarbons, it was shown how, when 
significant undercooling occurs on freezing, the observations can be 
analyzed to yield, in proper relation to the subsequent equilibrium 
part of the curve, the "zero" time, that is, the time at which crystalliza- 
tion would have begun in the absence of undercooling. In that paper, 
visual extrapolation of the equilibrium part of the curve, to "zero" 
time on freezing and to the corresponding time on melting, was used 
to determine the freezing point, defined as the temperature at which 

i This investigation is part of the work on Research Project 44 of the American Petroleum Institute from 
whose research fund financial assistance has been received. 

i Presented before the Division of Petroleum Chemistry at the meeting of the American Chemical Society 
at Cleveland, Ohio, April 5, 1944. 

3 Research Associate at the National Bureau of Standards, representing the American Petroleum Institute. 

* Figures in brackets indicate the literature references at the end of this paper. 
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an infinitesimal amount of crystals of the major component is in 
equilibrium with the liquid. 5 

The present paper describes a method, which may be used ana- 
lytically or graphically, for extrapolating the equilibrium part of a 
time-temperature freezing or melting curve to give the freezing point, 
T/, of the given hydrocarbon substance. There is also included a 
description of a method for deducing, within certain w r ide limits, from 
the same observations on the given substance if carried over a suffi- 
ciently large fraction of material crystallized or melted, the freezing 
point for zero impurity, T/ When the heat of fusion is known or can 
be estimated, the value of Tf — T f serves to give the amount of 
liquid-soluble, solid-insoluble impurity in the given substance. 

II. DERIVATION OF THE BASIC EQUATIONS 

The definitions, abbreviations, and symbols used in this paper will 
be exactly the same as those of the previous paper [1], except that the 
subscript / attached to the time 2, as in z f , denotes the time corre- 
sponding to the freezing point, T f . Some of the previous definitions 
and equations are repeated below for the sake of clarity and emphasis. 

For an ideal, or a sufficiently dilute solution, the thermodynamic 
relation between the temperature of equilibrium and the composition 
of the liquid phase, for the equilibrium between a liquid phase consist- 
ing of a major component and one or more other components and a 
crystalline phase of the major component alone, is [2, 3] 

-In iVx= -In (1-iVa) =^(1 + 1/2^2+ . . .) = 

(AHys^ . . .]. (i) 

In eq 1, the symbols are defined as follows: 

iVx = the mole fraction of the major component in the liquid phase; 
iV2=(l— iVi) = the sum of the mole fractions of all other com- 
ponents in the liquid phase; 
R=the gas constant, per mole; 

2^= the absolute temperature of the freezing point of the major 

component when pure (that is, when iVi = l or A^ 2 =0); 

AH/ =the heat of fusion, per mole, of the major component in the 

pure state at the temperature T fo ; 
ACp =the heat capacity per mole, of the pure liquid less that of the 
pure solid, for the major component in the pure state at 
the temperature T fo ; and 
T=the given temperature of equilibrium. 
The assumptions involved in the derivation of eq 1 are the following: 

(a) The liquid solution is ideal, or sufficiently dilute, so that the 
fugacity of the major component in the liquid phase is proportional 
to its mole fraction. 

(b) The heat of dilution is negligible over the range of concentration 
from Ni = 1 to the given value of N u 

5 It is important to note that this temperature, when properly measured, must have the same value whether 
approached from a higher temperature as in cooling or from a lower temperature as in warming. 
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(c) The heat of fusion of the major component in the pure state at 
the temperature T may be represented by a power series expansion, as 

&H f =AH f0 +AC Pe (T-T, )[l+ . . .]. (2) 

Introducing the following additional abbreviations, 

A=AHJRT, *, and B=llT,-ACJ2AH fo , (3) 

expanding —In (1— 7V 2 ), and rearranging, eq 1 becomes 

T f -T=^[\ + II2N 2 + . . .]-[B(T,-T)*+ . . .]. (4) 

Neglecting terms in iV 2 3 and higher, eq 4 becomes 

I7 -T=4 2 [1 + 1/22V 2 + . . .]_[z?(^ 2 ) 2 + . . .] (5) 

or 

T /0 -T=^[l + (l/2-B/A)N 2 + . . .]. (6) 

In a system of the type under investigation, in which the solid phase 
is restricted to the major component, let n±* represent the total num- 
ber of moles of the major component (in both the liquid and solid 
phases), n 2 represent the total number of moles of the other compo- 
nents (in the liquid phase), and n x represent the number of moles of 
the major component in the liquid phase under given conditions. 
Then the number of moles of the major component in the crystalline 
phase under the same conditions is %*— %. Let r denote the fraction 
crystallized of the total number of moles of all components, that is, 

r=(n 1 *-n 1 )l(n 1 *+n 2 ). (7) 

The mole fraction, in the liquid phase, of the sum of all components 
other than the major components is defined by 

N 2 =n 2 l(n x -\-n 2 ). (8) 

Let N 2 * represent the value of N 2 when no crystalline phase is present, 
that is, when 711=%*. Then 

iV 2 *=7i 2 /fal*+^2). (») 

Combination of eq 7, 8 , and 9 yields the following equation relating 
N 2 and r: 

N^NfKl-r). (10) 

Letting 

a=N 2 */A, and b=(l/2-B/A)(N 2 *) 2 /A } (11) 

and combining eq 10 and 6, there is obtained 

r^-rVuV ■ • ' = r »-rM'+^+ ■ 1 (12) 



■i-itwitf m-H urn* 
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It is possible to further simplify eq 12 by using the general relation 
that, when 



G^)"<>. 



a _ a r 6/a"]' 1 a r b/a 1 



(13) 



Assuming that the unwritten terms in eq 12 and 13 are negligible, or 
may be substituted one for the other without significant error, then 
eq 12 may be written 

In general, the fractional error introduced by making the approxi- 
mations leading to eq 14 will be of the order of 1/12 to 1/3 JV 2 2 , de- 
pending upon the value of B/A (see eq 3 and 4). That is to 'say, 
fractional errors in N 2f or correspondingly in T f —T, will be about 
1/3 to 1 1/3 percent when ^ = 0.2, and about 1/12 to 1/3 percent 
when N 2 =0.1. Equation 14 therefore would appear to be practically 
useful only when the amount of solute or impurity is below several 
tenths in mole fraction. 

The next step is to express r, the fraction crystallized or melted, 
in terms of the time. This is done by assuming that the experiment 
is one in which the rate of crystallization or melting of the major 
component is constant with time. 6 

Denoting the time by the symbol z, and letting z= z f when r=0, one 
has the relation 

r=k(z~z f ) t (15) 

where k is a constant characteristic of the experiment, positive in 
sign for a freezing experiment and negative for a melting experiment. 
The time at which freezing begins, or melting is complete, is z ff pro- 
vided thermodynamic equilibrium is maintained throughout the 
experiment. The determination of z f for an actual experiment has 
been discussed in reference [1]. 

Combination of eq 14 and 15 yields the following relation between 
the time and the temperature of equilibrium: 

T = T '-i-k*(z-z f y ( 16 ) 

where a' = a/(l — b/a) and k'=k/(l — b/a). 

Equation 16 relating T and z and eq 14 relating T and r are the 
equations of a rectangular hyperbola. The asymptote obtained by 
setting z= — <*> in eq 16 is T=T /q . 

Since there are three constants in eq 16, it is clear that three points 
on a portion of the curve satisfying eq 16 are sufficient to determine 
the values of the three constants, T /qJ a', and k'. However, if it is 

6 This is essentially the case in experiments performed with the apparatus described in reference fl] when 
the range of freezing (or melting) temperature is sufficiently small in comparison with the head of 
temperature. 
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desired to test the accuracy with which the data satisfies the equation, 
a greater number of points are necessary. 

The following theorem concerning a rectangular hyperbola will be 
required. Consider any four points (2\, z x ) } (T 2 , z 2 ), (T z , z 3 ) and 
T i} 2 4 ), on the hyperbola defined by eq 16. There are three ways 
in which two chords can connect these points, each point being the 
terminus of only one chord. These are illustrated in figure 1, in 
which m ]2 denotes the slope of the chord connecting points (T u z x ) 
and (T 2 , z 2 ), and similarly for ra 13 , ra 2 3, and m 14 . It is to be shown that 
the products of the slopes of the two chords is the same in each of 
the three cases 

mi 2 m u =m u m^=m u rn 23 . (17) 

Using eq 16, 



(T 2 -T x \ 



-a'V 



-[l-k'{z x -z f )][l-k'(z 2 -z f )] 



and 



m 12 m 34 . 



(a'k') 2 



■[l-k\z x -z f )][l-k\z 2 -z r )][l-k\z,-z f )][{l~¥{z,~z f )] 



(18) 



(18a) 



Equation 17 follows immediately, since the expression on the right- 
hand side of eq 18a is unchanged in value by any interchange of the 
subscripts 1, 2, 3, and 4. 




Time - 



Figure 1. — Construction to show three ways of connecting four points on a curve 
with two chords, with each point being the terminus of only one chord. 



III. DETERMINATION OF THE FREEZING POINT 

1. ANALYTICAL 

The procedure for determining the freezing point of the given 
substance analytically from the equilibrium portion of the time- 
temperature curve is illustrated by means of figure 2. Points G> H, 
and / are on the equilibrium part of the time-temperature curve, as 
defined by eq 18. The zero tune, at which crystallization would have 
begun in the absence of undercooling, has been determined [1], It is 
required to determine on the ordinate through the zrro time, ^ / , the 
freezing point T f . 
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Time — ■> 

Figuke 2. — Construction to show the analytical determination of the freezing point. 
From eq 17, it follows that 

/^ys^wfcsys^sx (19) 

\Zg—Zf/\Zi — Z h J KZi—Z/jyZn—Zg/ J 

Rearranging and using the relation 



one obtains 



\irV \lg—l hJ\Zi—z h J\z g —Zf/ 
Denoting by % t?, and w the dimensionkss ratios 

\z g -z f J 
there is obtained, on substitution into eq 21, the relation 

Equation 25 gives the value of T f , when there are known the tempera- 
ture and time corresponding to three points (as G, H, and I) on the 
equilibrium curve (as defined by eq 14) and the "zero" time, the time 



(20) 

(21) 

(22) 
(23) 
(24) 

(25) 
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at which crystallization would have begun in the absence of under- 
cooling. 7 

It is nearly always possible to select the point H equidistant in time 
between G and /, so that 0=1, Furthermore, it will usually be pos- 
sible to space the points 6 r , H, and / so that w is a simple ratio. Thus, 
if all three time intervals are equal 



(z g -z f ) = (z h -z g ) = (z i -z h y 

v=l 
w=3 

T f =T g +{T g -T x )l{Zu-\)) 



(26) 



When the first time interval is twice as great as the other two intervals, 
then 



(z g —Zf)=2(z h —z g )=2(z i —z h y 

v=l 

w=2 

T f =T g +(T g -T t )/(2u-l) 



(27) 



It is obvious that a large value of w corresponds to a short extrapola- 
tion to T f , and vice versa. 

Although the curve in figure 2 has been drawn to simulate a freezing 
curve, the above equations apply equally well to a melting curve, since, 
for the same substance, the melting curve is the mirror image of the 
freezing curve when the time rate of melting is made the same as 
the time rate of crystallization. 

2. GEOMETRICAL 

The manner of determining geometrically the freezing point, T f , 
is shown in figure 3. 

There are given the "zero" time, z f , and the time and temperature 
corresponding to the three points, G, H, and /, on the equilibrium 
part of the time- temperature curve defined by eq 16. It is required 
to locate the freezing point, T f , on the ordinate at the time z f . 

The s oluti on is as follows: 

Draw AB through I parallel to the time axis. 

Draw AC perpendicular to AB, at the time z f . 

Draw a line through points G and H, intersecting AB at E and 
ACatD. 

Draw a lin e through poin ts I and H, intersecting AG at J. 

Draw JK parallel to DE. 

Draw a line through points K and <?, intersecting AC at F. 

F is the desired point, representing the freezing point, T f . 

7 It is worthwhile remarking that eq 25 is not limited to the determination of the freezing point, Tf. 
Zf in eq. 24 may be any given time, and eq. 25 then yields the corresponding temperature, Tf. It has been 
implied in the foregoing discussion that the cooling curve before crystallization begins (or the warming 
curve after melting is complete) is sufficiently steep so that 2/ may be estimated exactly enough from a 
preliminary visual extrapolation to Tf. It is clear that if this is not the case, T/may still be obtained with 
any required accuracy by successive applications of eq 25, 



m m wwutmm-n 
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The proof, illustrated in figure 4, is as follows: 
Draw FL 

Draw NP through H, parallel to Aft. 
Then 

Jm/NP^NDINJ (28) 

But 

^U'lAK^NH/NT. (29) 

Therefore, 

ND/NJ==AI/AK. (30) 




H~ Time — * 

Figure 3. — Geometrical construction for determining T f and T/ Q . 
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But, by oq 17, 



(ND/NH)(AFfAI) = (AFlAK) (NJ/NH), 



which o'ivos 



in agreement with eq 30. 



ND/NJ=AIIAK, 



(31) 

(32) 



IV. DETERMINATION OF THE FREEZING POINT 
FOR ZERO IMPURITY 

The method described here for determining the freezing point for 
zero impurity, T f , can be properly applied only to those experiments 




Zf Time — * 

Figure 4. — Proof of geometrical construction for determining Tf and Tf Qt 
582610—44 4 
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in which the time-temperature observations extend over a sufficiently 
large fraction (probably of the order of }{ to %) of material crystallized. 
Furthermore, since that portion of the curve which is required for 
determining T fo is also the portion in which departures will most likely 
occur in the relations upon which eq 16 is based, it is seen that the 
determination of even an approximate value of T, will frequently be 
impossible for man}^ experiments. However, for completeness, the 
description of the procedure for evaluating T f by this method is given 
here so that it may be available for application to those experiments 
which do conform to the requirements. 

1. ANALYTICAL 

An equation for the freezing point for zero impurity, T f , may be 
arrived at as follows. Returning to eq 16, it will be noted that T, 
as defined there, approaches T f as z approaches — °° . This may be 
regarded simply as an analytical property of eq 16, so far as the present 
argument is concerned. 8 

Since eq 25 applies to any point on the curve defined by eq 16, we 
may imagine z f and T f in eq 24 and 25 to be replaced by z and T, 
respectively, corresponding to any point on the curve. Then as z 
approaches — oo f w approaches unity, and eq 25 becomes 

T,=T g +(^=p)- (33) 

o * \uv—l J 

Since w has been eliminated from eq 33, 2V does not depend on the 
location of z f . As pointed out following eq 25, it is nearly always 
possible to select the point H equidistant in time between G and I, 
so that 0=1. 

The precision of the determination of the freezing point for zero 
impurity rests heavily upon the value of the product (uv), which can 
be seen from eq 22 and 23 to be the ratio of the slope of the chord 
HI to the slope of the chord GH. 

2. GEOMETRICAL 

The manner of determining geometrically the freezing point for 
zero impurity, T f , is shown in figure 3. 

After having determined the freezing point, T f , as indicated in the 
preceding section, the procedure is as follows: 9 

Diaw a line through G parallel to AB and a line through / parallel 
to AC ) the two lines intersecting at L. 

Draw a line through points K and L, intersecting AC at M. 

M is the desired point, representing the freezing point for zero- 
impurity, T Tq . 

8 The extension of the time in eq 16 to -co may be easily pictured by considering a corresponding melting 
curve, which may be prolonged indefinitely in time by the continued addition of the pure major component 
in the solid phase to the existing solid phase as melting proceeds. This procedure serves to produce a liquid 
phase in which the number of moles of the major component approaches infinity, the composition of the 
liquid phase approaches Ni-l or iV*=0, and the temperature approaches the freezing point for zero im- 
purity, 7/ . 

The method of obtaining 7/ described here has the advantage that T/ 9 is obtained on the same vertical 
axis as Tf, so that the initial freezing point lowering (T/ — Tf) is immediately apparent. However, it 
should be pointed out that the value of TV obtained depends only on the points G, H, and I, and not at 
all on the time chosen for if. 
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The proof is as follows: 
Draw HI. 

Since point F is on the curve, eq 22, 23, and 33, for T f0 , may be 
rewritten with the three points G, H, and / replaced by F, G, and /. 

T *~ T M&=i} (34) 

Ht^t) (35) 



where 



\Zi-zJ 



(36) 



Equation 34 may be put in the form 

ffi>~!*Wv. (37) 



But from figure 4, 



\T f0 -T f J- 



AM/FM==IL/RL. (38) 



Also, from eq 35 and 36 and figure 4, 

u'v'=m OI /m FGy (39) 

where m G/ and m FG are the slopes of the chords ftFand FG, respectively. 
Combining eq 37, 38, and 39, 

IL/m GI =RLIm FG > (40) 

But both sides of eq 40 are equal to GL, from figure 4. 

V. UNCERTAINTY OF THE DETERMINATION OF THE 
FREEZING POINT AND OF THE FREEZING POINT FOR 
ZERO IMPURITY 

An expression may be derived for the error of T f or T fo , as calculated 
from eq 25 and 33, for any given experiment. For simplicity, one 
may assume that there is no error in z f , z gi z h , and z i} but that a 
uniform random error, a, is to be attached to T gy T h , and T t (fig. 2). 
The error in T f , oy, may then be calculated from the usual propaga- 
tion of error equation 

Equation 41 assumes that the errors in T^, T h , and T t are, on the 
average, uncorrelated, and gives the error for T f in terms of the 
average of these errors, a. 

The derivatives in eq 41 may be obtained by differentiation of 
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eq 25, noting that u, but not v and w, is a function of T Wi T h , and T t . 
They may be put in the following form (with the aid of eq 22): 

(bTA u*vw(vw+l) .. 

\w g / (uvw-i) 2 ~ VU) 

fdT f \_ vwiu+iy 

\bTJ~ (wyw--\y { * 6) 

\bTj {uw)-\) 2 K%v 

When eq 42, 43, and 44 are substituted into eq 41, there is obtained 
the relation: 

(a A 2 _ (vw+ l)\u A v 2 w 2 + l)+v 2 w*(u+ 1) 4 

\a) (uvw-iy (4 ' J) 

which expresses (a/fa) as a function of u, v, and w. Equation 45 
may be put in a more useful form as follows: 

For any given segment of an equilibrium curve, the greatest 
accuracy m T f (minimum oy) will be obtained by choosing z and z t 
at the beginning and end of the segment in question (corresponding 
to a maximum value of w). With z § and z t fixed, the following 
relation then holds: 

uv=c, (46) 

where c is a constant characteristic of the given segment of the curve, 
and independent of the choice of z h . Equation 46 may be derived in 
the following way. 

It is evident from eq 22 and 23 that the product (uv) is t he ratio 
of the slope of the chord /// to the slope of the chord OH (fig. 2) 
From eq 18 c is then given by 



' \l-k%)' 



(47) 



which is independent of z h . 

If one eliminates u from eq 45 and 46, there is obtained the relation 



(?)■= 



v 2 (cw-iy 

A numerical investigation of eq 48 shows, as would be expected, 
that, for the range of conditions encountered in practice, the optimum 
value of v is in the neighborhood of unity (corresponding to z h equi- 
distant between z g and z<). When v=l, u=c, and eq 48 becomes 

fa f \ 2 _ (w+ 1) 2 (^W+ 1) + w 2 (u+ iy _ 

W (uw-iy ^ } 

The error in T f , oy , may be obtained immediately by setting w 
equal to unity in eq 48 or 49. When v~l, 

_4(uH-i) + M-i) 4 



(vT 



(«-!)• (60 » 



mm 
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From eq 49, (07/0-), for v=l, has been calculated for w=l, 1.5, 2, 
3, 5, and co ; and for u—l } 1.25, 1.5, 2, 3, 4, and ». The results are 
presented in figure 5, in which the abscissas represent (lfu) 9 and the 
ordinates represent (V//<r), on a logarithmic scale. The curve for 
w=l represents {a f /a) } and rises to infinity at (l/w)=l, which cor- 
responds to a freezing curve of constant slope. The other curves, 
representing (07/0-) for various values of w } are finite for all values 
of (1/u). 
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Figuke 5. — Chart to show the uncertainty of the determination of T/ and T/ Q . 

Figure 5 is expected to be convenient for estimating the uncertain- 
ties of T f and T f in practice. In calculating T f and T f from eq 25 

and 33, u and w are first evaluated from eq 22 and 24 (v=l). The 
values of (07/ o-) and (<r//(r)may then be read immediately from figure 5. 

It remains for the experimenter to estimate a reasonable value for <r. 
In any case, figure 5 indicates relative errors for various values of 
u and w. 

The errors discussed above do not include those arising from 
departures from the ideal solution laws or from lack of constancy of 
the rate of crystallization or melting (see appendix, page 212.) 
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VI. ALTERNATIVE METHOD FOR DETERMINING THE 
FREEZING POINT (AND THE FREEZING POINT FOR 
ZERO IMPURITY) 

The foregoing methods for determining T f (and, where appropriate, 
Tf ) may be used when the experimenter is able from previous expe- 
rience to select with some assurance the equilibrium portion of the 
time-temperature freezing or melting curve. When there is some 
doubt as to what portion, if any, of the time-temperature freezing or 
melting curve constitutes the equilibrium portion, the following 
procedure may be used to select that portion of the curve that is 
consistent with the hyperbolic form of eq 16. At the same time, 
there is obtained a value for the freezing point, T f , and within certain 
wide limits, a value for the freezing point for zero impurity, T/ , if 
the observations cover a sufficiently large fraction of the material 
crystallized or melted. 

Let (T m , 2 m ) be an arbitrary fixed point chosen on the equilibrium 
portion of a time- temperature curve. Then it may be shown from 
eq 16 that a plot of each of the following functions against the varia- 
ble indicated will yield a straight line: 



(^) versu <^;) 

versus (T-T m ), or T 



( T-T m \ 

\z—z m ) 

( Z — Zm \ 

\T-TJ 



versus (z—z m ), or z 

The plot of (z—z m )/(T—T m ) versus (z— z m ) is the most convenient 
for the purposes of this paper. Define new quantities y and x as 
follows: 

x=(z—z m ); x f =(z f —z m ) (52) 

Then it follows from eq 16 that 

y=Vo+sx; y f =y +sxf (53) 



where 



(5tf 



s -(t, -tJ- 

Thus a plot of y against x yields a straight line for that portion of 
the curve, presumably the equilibrium portion, which is consistent 
with eq 16. Extrapolation of the straight line to x=x f (z=Zf), 
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yields y f ; T f may then be calculated from the equation 



similarly, 



T,= T m +% (56) 

Vt 

T f ^T m +l/s=T m +^^, (57) 

Xf 

since x m =0.^ The smoothed temperature at any other time is given 
by the relation 

T=T f -^+~- (58) 

The point (T m , z m ) should be selected near the middle of that por- 
tion of the time-temperature curve which the experimenter regards 
as most likely to represent equilibrium. The linearity of the y versus 
x plot is then an indication of the range on either side of (T m , z m ) in 
which equilibrium has apparently been maintained. It should be 
noted that the error in the calculated value of y increases rapidly as 
x decreases. In plotting y it is better to plot circles, of radii inversely 
proportional to x, than points. This will insure approximately the 
correct weighting of the points in drawing the straight line. It 
should also be noted that a slight displacement of ( T m , z m ) from the 
true curve will have a negligible effect on the location of the straight 
line, provided (T m , z m ) is about in the middle of the equilibrium por- 
tion of the curve. 

Equation 57 relates T f0 to the slope, s, of the y versus x plot. How- 
ever, it is clear that s is a sort of second derivative of the time-tem- 
perature curve. It is to be expected, therefore, that the value of T fQ 
calculated from eq 57 may be quite inaccurate. It is possible to re- 
verse the procedure, if T fo is known, and calculate the value of the 
slope of the y versus x plot to be expected on the basis of eq 16. 

The y versus x plot is principally of use in extrapolating to obtain 
T f} by eq 56, and also in obtaining smoothed values of the tempera- 
ture at any other time, by eq 58. The smoothed values of the tem- 
perature may be used in conjunction with the other methods of 
obtaining T fQ as described [1]. 

[1]. 

VII. DISCUSSION 

In connection with the practical application of the foregoing 
methods for determining from time-temperature freezing and melting 
curves the freezing point, T f , and the freezing point for zero impurity, 
T f0 j the following points are important: 

(a) The experimenter must assure himself that the time-tempera- 
ture relations substantially follow eq 16 over that part of the experi- 
ment which represents thermodynamic equilibrium between the 
liquid and solid phases. (See section VI and the appendix of this 
paper.) 

(b) The precision of determining the freezing point, T f , will be 
nearly as good as the determination of temperatures along the equi- 
librium part of the curve, except when the time and temperature 
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covered by undercooling are very large. In terms of the uncertainty 
in the temperatures along the equilibrium part of the curve, the un- 
tainty in T f is given by eq 49 and illustrated in figure 5. 

(c) The precision of determining the freezing point for zero impur- 
ity, T f , by the method described, will be much less precise than the 
determination of the freezing point. In terms of the uncertainty in 
the temperatures along the equilibrium, portion of the curve, the un- 
certainty in T f() is given by eq 50 and illustrated in figure 5 by the 
curve for which w=l. The practical utility of determining T f() by 
this method depends entirely upon the equilibrium portion of the 
curve covering a sufficiently large range of r (fraction crystallized, or 
melted) so that the product uv in eq 33, the ratio of the slopes of the 
two chords, will have a value appreciably greater than unity. If the 
observations of the given experiment do not cover a sufficiently large 
fraction of material crystallized or melted, no useful determination of 
T fQ may be made by this method. For any given experimental pro- 
cedure, it is desirable to check the value of T fQ obtained by this method 
with the value obtained by the method previously described [1], which 
has been shown to give reliable results. 
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IX. APPENDIX 

In the foregoing discussions, it has been assumed for simplicity that the rate of 
crystallization is a constant, or that r is linear in z, as given by eq 15. It is of 
interest to consider the effect of a varying rate of crystallization. Two types of 
variation may be distinguished: 

(a) The rate of crystallization varies initially, but approaches a constant value 
as the experiment proceeds. 

(b) The rate of crystallization varies progressively as the experiment proceeds, 
and does not tend to approach a constant value. Supercooling should give rise 
to a variation of type (a). Variations of type (a) will destroy, for a certain time, 
the rectangular hyperbolic relation between Tand z, eq 16, on which the linearity 
of the {z—Zrn)l{T-T m ) versus(z— z m ) plot depends. It may therefore be assumed 
that in the linear region of this plot such variations have decreased sufficiently to 
be neglected. 

Variations of type (b) are of more importance, since they increase throughout 
the experiment. Such variations would be caused, for example, by any factor, 
such as a decreasing (or increasing) head of temperature, that alters the rate of 
heat transfer and hence the rate of crystallization (or melting). 

If the variation is of type (b), r may be represented by a power series in {z — z f ): 

r==k ( z - Zf )[l- e ( z - Zf ) + . . .]. (59) 

Now consider the following expansion, valid for [e(z — z/)] 2 less than unity: 

[l^ll^l = k { z~z f )[l-e(z-z f )+ . . .]. (60) 

Ll + e(z — zf) J 

Therefore, when the higher order terms in eq 59 and 60 may be neglected, r may be 
written: 

Mz — z f ) 
r =l + e(z-z f )' < 61 > 
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Equation 61 is a rectangular hyperbolic relation between r and z. Furthermore, 
eq 14 is a rectangular hyperbolic relation between T and r. It is easily shown 
that the relation between Tand z obtained by substitution of eq 61 in eq 14 is also 
rectangular hyperbolic. Specifically, eq 16 is replaced by the equation: 

T-(t i a'e \ a'k'/(k'-e) ( . 

T -V f «+k'-e) l-(k'-e){z-z f ) (62) 

Equation 62 reduces to eq 16 for e equal to zero. 

It is clear from eq 62 that all the results previously derived on the basis of eq 16 
are unchanged, provided k' is replaced by (k r — e), a' is replaced by [a'k '/ (k' — e)], 
and T/ is replaced by [Tf Q -\-a'el(k r -e)}. Of the previous results which are unaffec- 
ted, there may be mentioned the theorem concerning the slopes of the chords, eq 
17, and the linearity of the {z—z m )/(T—T m ) versus (z— z m ) plot. The calculated 
value of Tf will not be in error. However, the calculated value of 7> will be in 
error by the amount [a'el(k f — e)], corresponding to the fact that the asymptote of 
the hyperbola is shifted. 

Washington, January 6, 1944. 
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